Event-by-event simulation of a quantum eraser experiment* 
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We present a computer simulation model that is a one-to-one copy of a quantum eraser experiment 
with photons (P. D. D. Schwindt et ah, Phys. Rev. A 60, 4285 (1999)). The model is solely based 
on experimental facts, satisfies Einstein's criterion of local causality and does not require knowledge 
of the solution of a wave equation. Nevertheless, the simulation model reproduces the averages 
as obtained from the wave mechanical description of the quantum eraser experiment, proving that 
it is possible to give a particle-only description of quantum eraser experiments with photons. We 
demonstrate that although the visibility can be used as a measure for the interference, it cannot be 
used to quantify the wave character of a photon. The classical particle-like simulation model renders 
the concept of wave-particle duality, used to explain the outcome of the quantum eraser experiment 
with photons, superfluous. 
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I. INTRODUCTION 

According to wave-particle duality, a concept of quan- 
tum theory (QT), photons exhibit both wave and par- 
ticle behavior depending upon the circumstances of the 
experiment 1 ". The wave and particle behavior of pho- 
tons is believed to be complementary. When we know 
(observe) the which-way (WW) information (particle be- 
havior), there is no interference pattern (wave behavior) 2 . 
Parameters quantifying the interference and the WW in- 
formation are the visibility V and the path distinguisha- 
bility T>, respectively. According to the complementarity 
relation of QT, V 2 + V 2 < l&i 

In 1982, Scully and Driihl proposed a photon interfer- 
ence experiment, called "quantum eraser"—, in which the 
photons are labelled by WW markers (three-level atoms). 
In this experiment, we know (but not observe) the WW 
information of the photons and then we expect that there 
is no interference. However by erasing the WW infor- 
mation afterwards by a "quantum eraser" , the interfer- 
ence pattern can be recovered^. The interference pattern 
can even be recovered after the data have already been 
recorded and saved in a file£. 

Quantum eraser experiments have been described "as 
one of the most intriguing effects in quantum mechanics" , 
but have also been regarded as "the fallacy of delayed 
choice and quantum eraser"—. Clearly, they challenge 
the point of view that the wave and particle behavior of 
photons are complementary: The observation of interfer- 
ence, commonly associated with wave behavior, depends 
on the way the data is analyzed after the photons have 



FIG. 1: Schematic diagram of the experimental setup for the 
quantum eraser experiment with photons studied in Ref. H|. 
BS: beam splitter; PBS: polarizing beam splitter; HWP0 and 
HWP1: half- wave plates; QWP: quarter- wave plate; Do, D\: 
detectors; ci: phase shift introduced in Pathl. 



passed through the interferometer. 

The question that we answer in the affirmative in this 
paper is: "Can we simulate a quantum eraser experiment 
without invoking concepts of quantum theory and with- 
out first solving the wave mechanical problem?" 



A. Quantum eraser with photons 

1. Experimental realization 



The quantum 
several different 
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eraser has been implemented in 
experiments with photons, atoms, 
etc-ii— — — — — — . Although much more difficult to real- 
ize experimentally, quantum erasers may also be realized 
with quantum dot o 14 ' 15 and mesoscopic electromechani- 
cal devices 1 ^. 

In Ref. 8 , Schwindt et al. reported an experimental re- 
alization of a quantum eraser in which the polarization 
of the photons has been used to encode the WW infor- 
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mation. In this paper, we focus on this particular exper- 
iment. The experimental setup (see Fig. [T]) consists of 
a linearly polarized single-photon source (not shown), a 
Mach-Zehnder interferometer (MZI) of which the length 
of Pathl (see Fig. [T]) can be varied, inducing a relative 
phase shift <fi between PathO and Pathl, an adjustable 
analysis system which is a combination of a quarter- 
wave plate (QWP), a half-wave plate (HWP) HWP1, and 
a calcite prism operating as a polarizing beam splitter 
(PBS). Another adjustable HWP, HWPO, is inserted in 
PathO of the MZI to entangle the photon's path with its 
polarization. 

According to Rcf. 8 the pictorial description of the ex- 
perimental observations is as follows. If a photon, de- 
scribed by a pure, vertically polarized state V is injected 
into the interferometer with the HWPO set to 45°, then 
the photon that arrives at the second beam splitter (BS) 
of the MZI carries a WW marker: The photon is in the 
horizontally polarized state H if it followed PathO and 
it is in the V state if it followed Pathl. If the opti- 
cal angle of HWP1 is zero, there will be no interference 
(V = 0) and the detectors give us the full WW infor- 
mation of each detected photon (T> = 1). If the optical 
angle of HWP1 is nonzero, the H and V states interfere 
(0 < V < 1) and the WW information of each photon 
will be partially or completely "erased" (0 < T> < 1). 
Thus, by varying the optical angle 6hwp\ oi HWP1, the 
illusion is created that the character of the photon in the 
MZI "changes" from particle to wave and vice versa. If 
photons described by a completely mixed, that is an un- 
polarized, state are emitted, then no WW information 
can be obtained and also no interference can be observed 
(T> = V = 0), independent of the orientation of HWPO. 
However, varying 9hwpi can still lead to a recovery of 
interference (0 < V < 1). For photons described by a 
partially mixed state, a state that can we expressed as 
containing a completely mixed component and a pure 
component, partial WW information can be obtained. 
Since the completely mixed component contains no WW 
information and displays no interference, the maxima of 
T> and V are smaller than one and numerical equal to the 
state purity. Also in this case complete visibility can be 
recovered by varying 9hwpi ■ 



2. Event-by-event simulation model 

It is important to realize that the counter-intuitive 
features of quantum eraser experiments result from at- 
tempts to apply the concepts and the formalism of QT 
to a description of the experimental results in terms of 
individual events^. Logically speaking, there are two pos- 
sibilities: 

1. We accept the postulate that it is fundamentally 
impossible to give a logically consistent description 
of the experimental results in terms of individual 
events, that is we accept that there is no expla- 



nation that goes beyond the quantum theoretical 
description in terms of averages over many events. 



2. We search for an explanation of the experimental 
facts that goes beyond a description in terms of 
averages. 



In this paper, we demonstrate that the second option is 
a viable one. Thus, we adopt the point of view that al- 
though QT correctly predicts averages of many detection 
events, it has nothing to say about individual events^. 

We propose an event-by-event simulation model that 
is a one-to-one copy of the quantum eraser experiment 
reported in Ref. 8 . The simulation model describes a 
particle-like, classical, local and causal dynamical sys- 
tem. Each component of the laboratory experiment such 
as the single-photon source, the BS, HWP, QWP, and 
PBS are simulated by corresponding algorithms. By con- 
necting the output (s) of one component to the input (s) 
of another one, we construct the simulation equivalent 
of the experimental setup depicted in Fig. [TJ By con- 
struction this network of dynamical systems satisfies Ein- 
stein's criterion of local causality. The data is analyzed 
by counting the detection events, just as in the real ex- 
periment. 

We demonstrate that our model reproduces the re- 
sults of QT, that is the averages predicted by QT and 
confirmed by experiment 8 , without first solving a wave 
equation. In fact, we show that it is possible to give an 
entirely classical, particle-only description for the single- 
photon quantum eraser experiment reported in Ref. 8 . We 
show that the interference patterns, commonly associated 
with wave behavior, can be built up by many particles 
having full WW information (we can always track the 
photons during the simulation) that arrive one-by-one at 
a detector. 

The work of this paper builds on earlier 

work 17^ i 19 1 20 i 21 1 22 i 23^4 1 25^6 1 27^8 1 29 ^ dem0 nstrateS 

that quantum phenomena can be simulated on the 
level of individual events without first solving a wave 
equation and even without invoking concepts of QT, 
wave theory or probability theory. Specifically, in 
our earlier work we have demonstrated that it is 
possible to simulate event-by-event, a single-photon 
beam splitter and Mach-Zehnder interferometer exper- 
iments, Einstein-Podolsky-Rosen-Bohm experiments 
with photons, Wheeler's delayed choice experiment 
with single photons, the double-slit and two-beam 
interference, quantum cryptography protocols, and 
universal quantum computation. The latter proves 
that in principle we can perform an event-by-event 
(particle- like) simulation of any quantum system 30 . 
Some interactive demonstration programs are available 
for downloa d 31 ' 32 ' 33 . 
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B. Irrelevance of Bell's theorem 

It is not uncommon to find in the recent liter- 
ature, statements that it is impossible to simulate 
quantum phenomena by classical processes. Such 
statements are thought to be a direct consequence of 
Bell's theorem 3 ^ but are in conflict with other work 
that has pointed out the irrelevance of Bell's theorem 

35 ,36.37,38.39,40.41,42,43,44,45,46,47,48,49.50,51.52.53,54,55,56.57 

A survey of the literature suggests that, roughly speak- 
ing, physicists can be classified as those who believe 
in the reasonableness of Bell's arguments, those who 
advance logical and mathematical arguments to show 
that a violation of Bell's (and related) inequalities does 
not support the far-reaching conclusions of the former 
group of physicists and those who do not care about 
Bell's theorem at all. The authors of this article belong 
to the second group. 

Although we expect discussions of philosophical or 
metaphysical aspects of Bell's theorem to continue for- 
ever, as explained in a review article that has appeared in 
this journal 2 —, from the viewpoint of simulating quantum 
phenomena on a digital computer, Bell's no-go theorem 
is of no relevance whatsoever. 

This conclusion is supported by several explicit exam- 
ples that prove that it is possible to construct algorithms 
that satisfy Einstein's criteria for locality and causality, 
yet reproduce exactly the two-particle correlations of a 
quantum system in the singlet state, without invoking 
any concept of quantum theory 2 ^ 2 ^ 3 -! 2 ^ 2 ^^. It is there- 
fore an established fact that purely classical processes 
can produce the correlations that are characteristic for a 
quantum system in an entangled state, thereby disposing 
of the mysticism that is created by Bell's no-go theorem. 

The key point is to realize that QT or the probabilistic 
models proposed by Bell cannot, on a fundamental level, 
address the (non)existence of algorithms, that is of well- 
defined processes, that give rise to the distributions of 
the events, described by these theories/models. 

The philosophy behind our simulation approach is very 
simple: If we can construct an algorithm that 

1. does not rely on the solution of a wave equation, 

2. satisfies the elementary criteria of locality and 
causality as formulated by Einstein, 

3. produces data of the same type as the data col- 
lected in the laboratory experiment, 

4. by analyzing the simulated data according to the 
procedure used to analyze the experimental data 
leads to the same conclusion, namely that certain 
averages of the raw data agree with the quantum 
theoretical description of the whole experiment, 

5. contains algorithms that simulate the various com- 
ponents (beam splitter, etc.) of the experiment and 
can, with no change, be re-used to simulate other 
experiments, 



then we may conclude that we have built a simulation 
model for the laboratory experiment. 

Loosely speaking, if the experimenter would be unable 
to distinguish between data recorded in a genuine exper- 
iment and data provided by the simulation algorithm, 
then the experiment has been "de-mystified" in the sense 
that we have found a process that offers a description of 
the observed phenomena on the level of individual events 
and without invoking (concepts of) wave theory. 

To avoid possible misunderstandings, the work pre- 
sented in this paper is not concerned with an interpreta- 
tion or an extension of QT nor does it affect the validity 
of QT as such. QT describes the collective result of many 
events, that is averages of many events, extremely well 
but does not provide a description on the level of indi- 
vidual clicks of a detector—. 



C. Structure of the paper 

Section [II] reviews the standard concepts of QT that 
are needed to give a quantum theoretical treatment of 
the quantum eraser experiment 8 . Section lllll discusses 
the general ideas that underpin our event-by-event simu- 
lation approach. We address the fundamental problem of 
reconciling the observation of "clicks" with a wave me- 
chanical theory from the viewpoint of algorithms, pro- 
cesses and computation. We show that in general, it is 
impossible to attribute "clicks" to individual wave ampli- 
tudes and explain how our simulation approach circum- 
vents this fundamental problem. Section [TV1 explains how 
the pure and mixed states of a quantum systems can be 
represented in our simulation approach. In SectionfVl we 
specify the simulation model in full detail. Data of event- 
by-event simulations of the quantum eraser experiment 
are presented in Section EH We show that our classical, 
particle-like simulation model reproduces all the results 
of QT for this experiment. Our conclusions can be found 
in Section lylTl 



II. QUANTUM THEORY 

In QT, a system is described by the state \ct), a vector 
in a Hilbert spaced. This vector can be written as a lin- 
ear combination of a complete set of orthonormal basis 
states \i) for i = 1, . . . , d where d denotes the dimension 
of the Hilbert space. These basis states are chosen such 
that they facilitate the formulation of the model. The 
amplitude for a quantum system to go from a state \a) 

to another state is given by (j3\a) = Ylt=i (/^K) (*l a )- 
With respect to the basis states {\i}}, the optical appa- 
ratus T is defined through its transition matrix elements 
(i\T\j). If the optical apparatus T induces a transition 
from the state |a) to the state the amplitude for this 
transition is given by {P\T\a) = J2i,j=i O 3 !*) {A T \j) 
Finally, the probability Prob(/3, a) for this transition to 
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occur is related to the amplitude through the Born rule 
Prob(/3, a ) = |(/3|7»| 2 . (1) 

According to the above scheme, we can easily calculate 
the predictions of QT for the experiment shown in Fig.[T] 
The basis states correspond to H or V polarized photons 
that travel along PathO or Pathl. The transition matrices 
of the optical components such as the BS, PBS, HWP 
and QWP can be found in Rcf. 60 < and in the appendix. 
In the appendix, we also give the quantum theoretical 
expressions for the visibility for the experiment depicted 
in Fig. [T] that will be used for the comparison with our 
simulation results. 

The above formulation assumes that the quantum sys- 
tem is in the pure stated. Some of the experiments re- 
ported in Ref.— require a description in terms of a mixed 
stated. A system is in a mixed state if it is in one of 
its m pure states |ck2), ■ ' ■ ,\ a m) with probability 
Pi,P2j""' ,Pm, respectively^. A quantum system in a 
mixed state is conveniently described through the den- 
sity matrix-^ 



ABA 



P = ^Pj\ ot i){° t i 
3=1 



(2) 



where it is assumed that Y^j=iPj = 1; Pj ^ f° r 3 — 
1, . . . , m, and that the states \aj) are normalized such 
that Trp = 1. According to QT, for a system in a mixed 
state p, the expectation value of the operator Q is given 



(n)=Tt P Sl = '£ i p j {a j \n\a j 



(3) 



3=1 



III. EVENT-BY-EVENT SIMULATION 

Our event-based simulation approach is unconven- 
tional in that it does not require knowledge of the wave 
amplitudes obtained by first solving the quantum the- 
oretical problem nor do we first calculate the quantum 
potential (which requires the solution of the Schrodinger 
equation) and then compute the Bohm trajectories of 
the particles. Instead, the detector clicks are generated 
event-by-event by locally causal, adaptive, classical dy- 
namical systems. Our approach employs algorithms, that 
is we define processes, that contain a detailed specifica- 
tion of each individual event which, as we now show, 
cannot be derived from a wave theory such as QT. 

To understand the subtleties that are involved, it is 
helpful to consider a simple example. Let us consider the 
MZI unit of the quantum eraser and omit the polarization 
label of the photons. According to QT, the amplitudes 
bo and b\ to observe a photon in PathO or Pathl after the 
second BS are related to the input amplitudes ao and a\ 
by£I 



bo 
bi 



1 i 
i 1 







1 i 
i 1 



o 

«1 



Of) 

«1 



(4) 



Let us assume that ao — I and ai = 0, meaning that the 
photons enter the MZI through PathO only. The proba- 
bilities Po (Pi) for a click in detector Dq (Pi) are given 
by 



Pk = 



53 X! A k,j B 3,i A i;0 



k = 0,1. 



(5) 



Using Eqs. ((3]) and (JSJl a simple calculation yields a 
closed form expression for P}.. Once we know Pk, it is 
trivial to use it as input for a process that generates clicks 
of the detectors Do and D\. This approach relies on 
what we call the "solution" of the quantum theoretical 
problem. It is irrelevant whether we have a closed form 
expression for Pk or only know Pk in tabulated form. 
The point is that we analytically worked out the sums 
over the indices i and j in Eq. ([5]). Let us now assume 
that we do not know how to perform the sums over the 
indices i and j in Eq. (J5J by ourselves and that there is 
some "magical process" that carries out the sum for us. 
In other words, we assume that we do not know Po and 
Pi. 

In practice, any process that performs the sums in 
Eq. ((51) by selecting (one-by-one) the pairs from 
the set S = (0,0), (1,0), (0, 1), (1, 1) defines a sequence of 
"events" The key question now is: Can we identify 

the selection of the pairs with "clicks" , events registered 
by a detector? We now prove that this is impossible. 

A characteristic feature of all wave phenomena is that 
not all contributions to the sums in Eq. ([5]) have the same 
sign: In wave theory, this feature is essential to account 
for destructive interference. But, at the same time this 
feature forbids the existence of a process of which the 
"events" can be identified with the clicks of the detector. 

This is easily seen by considering a situation for which, 
for instance, Po = 0. In this case, the detector Po should 
never click. However, according to Eq. l[5|). the process 
that samples from the set S produces "events" such that 
the sums over all these "events" vanishes. Therefore, if 
we want to identify these "events" with the clicks that we 
observe, we run into a logical contradiction: To perform 
the sums in Eq. ([5]), we have to generate events that 
in the end cannot be interpreted as clicks since in this 
particular case no detector clicks are observed. 

Thus, the conclusion is that the individual terms in 
expression Eq. §5§ do not contain the ingredients to define 
a process that generates the clicks of the detectors that 
we observe. 

The crux of our event-by-event simulation approach is 
that we do not start from expression Eq. |(5J) but con- 
struct a process that converges to Eq. ((SJ) while gen- 
erating events that correspond to the observed events. 
To grasp this idea, consider the well-known Metropo- 
lis Monte Carlo (MMC) method for solving statistical 
mechanical problem a 62 ' 63 . The MMC method generates 
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states S, events in our terminology, with a probability 
density^^ 

e -E(S)/k B T 
P ( S } = ^ s g-_E(S)/fc B T' ( 6 ) 

where E(S) denotes the energy of the state S, ks is 
Boltzmann's constant and T is the temperature. At first 
sight, sampling from Eq. ^ is impossible because in all 
but a few nontrivial cases for which the partition function 
J2s e~ E ^ s ^ kBT is known, we do not know the denomina- 
tor. MMC solves this problem by constructing a Markov 
chain that generates a sequence of events S such that 
asymptotically these events are distributed according to 
the (unknown) probability density Eq. ([5]) 62 iS. 

The analogy with our event-by-event simulation ap- 
proach is the following. Although very different in all 
technical details, our event-based method uses a deter- 
ministic process of which the sampling distribution con- 
verges to the unknown (by assumption) probability dis- 
tribution Pk for k = 0,1. Initially, the system does 
not know about this limiting probability distribution and 
hence, during a short transient period, the frequencies 
with which events are generated may not correspond to 
this distribution. However, for many events, which is the 
situation described by QT, these first few "wrong" events 
disappear in the statistical fluctuations and are therefore 
irrelevant for the comparison of our event-based simula- 
tion results with QT. It should be clear that the foregoing 
does not depend on the specific example that we used for 
the purpose of illustration. 

Let us now discuss the general aspects of our simula- 
tion approach. The simulation algorithms that we con- 
struct are most easily formulated in terms of events, mes- 
sages, and units that process these events and messages. 
Taking the quantum eraser experiment as an example, 
in a pictorial description, the photon is regarded as a 
messenger, carrying a message that represents its time- 
of-flight (phase) and polarization. In this pictorial de- 
scription, we may speak of "photons" generating the de- 
tection events. However, these so-called photons, as we 
will call them in the following, are elements of a model 
or theory for the real laboratory experiment only. The 
only experimental facts are the settings of the various 
apparatuses and the detection events. What happens in 
between activating the source and the registration of the 
detection events belongs to the domain of imagination. 

The processing units mimic the role of the optical com- 
ponents in the experiment and the network by connecting 
the processing units represents the complete experimen- 
tal setup. The standard processing units consist of an 
input stage, a transformation stage and an output stage. 
The input (output) stage may have several channels at 
(through) which messengers arrive (leave). Other pro- 
cessing units are simpler in the sense that the input stage 
is not necessary for the proper functioning of the device. 
A message is represented by a set of numbers, conven- 
tionally represented by a vector. As a messenger arrives 
at an input channel of a processing unit, the input stage 



updates its internal state, represented by a vector, and 
sends the message together with its internal state to the 
transformation stage that implements the operation of 
the particular device. Then, a new message is sent to the 
output stage which selects the output channel through 
which the messenger will leave the unit. At any given 
time, there is only one messenger being routed through 
the whole network. There is no direct communication be- 
tween the messengers. From this general description, it 
should already be clear that the process that is generated 
by the collective of classical dynamical systems is locally 
causal in Einstein's sense. Our simulation approach docs 
not rely on concepts of probability theory but instead, it 
generates events by way of classical, dynamical processes, 
the frequencies of events of which converge to the quan- 
tum theoretical results as the dynamical system relaxes 
to its stationary state. 



IV. SIMULATION OF PURE AND MIXED 
STATES 



In QT, the pure state is a description of the whole 
experiment, not of the individual events that are recorded 
by the detectors^^. In our simulation approach, the 
messages carried by the messengers represent the pure 
state, corresponding to a density matrix of the form p = 
\&k) ( a k\, that is pj = for all j ^ k and pf. = 1. In our 
simulation approach, the messages are constructed such 
that a large collection of them yields the same averages as 
those we obtain from quantum theory. Loosely speaking, 
we may say that a set of N (TV sufficiently large) messages 
of a certain type correspond to a pure state. 

In the more general case, QT describes the whole ex- 
periment through the mixed state Eq. ([2]). We simu- 
late the mixed state by the following procedure. Given 
Pi,... ,Pm> we pick an index k S {1, . . . , m} using a 
pseudo-random number and then send Nk messages of 
type k (corresponding to the pure state \oik)) through 
the network of processing units that represent the quan- 
tum system. The precise value of Nk is unimportant, as 
long at it is large enough to let the classical dynamical 
system mimic the pure state \ak}- 

For the case at hand, the quantum eraser, the source 
can emit a pure state, a linear combination of V and H 
polarized photons, or it can produce a mixed state of 
the two^. Thus, we have m = 2 and it what follows we 
will label the TV's by the subscripts V and H to facil- 
itate the comparison with the terminology used in the 
experiment^. Although not essential, in our simulation 
we simply choose Ny = Nu and denote the probabilities 
for the V- and iJ-polarized photons in a mixed state by 
Pv and ph, respectively. 



6 




cosy/, n 
sin y/[' n 

C0S Wi„ 
sin ic, 
f- 



DLM 


T 



















(c,s:.) 

n 













(M) 




T 
















1 




(CO 





















FIG. 2: Diagram of a DLM-based processing unit that per- 
forms an event-based simulation of a beam splitter (BS). 
The processing unit consists of three stages: An input stage 
(DLM), a transformation stage (T) and an output stage (O). 
The solid lines represent the input and output channels of 
the BS. The presence of a message is indicated by an arrow 
on the corresponding channel line. The dashed lines indicate 
the data flow within the BS. The transformation matrix T is 
given in Eq. (115[l , 



field components in the plane orthogonal to k§i. Adding 
another clock to the messenger suffices to model the sec- 
ond electric field component orthogonal to the first one, 
and hence the fully polarized wave^l. 

Thus, each messenger carries a message represented by 
a six-dimensional unit vector 



Yk,l 



( C0S Ki \ 

COs£fc,i 

V sm£ fc ,z ) 



(7) 



where the superscript H (V) refers to the horizontal (ver- 
tical) component of the polarization and ip^i, 4>Xii an d 
£,k,i represent the time of flight and polarization of the 
photon, respectively. It is evident that the representation 
used here maps one-to-one to the plane- wave description 
of a classical electromagnetic field^, except that we as- 
sign these properties to each individual messenger, not 
to a wave. The subscript I > numbers the consecutive 
messages and k = 0, 1 labels the channel of the BS at 
which the message arrives (see below). 



V. SIMULATION MODEL 

As explained earlier, our simulation algorithm can be 
viewed as a message-processing and message-passing pro- 
cess: It routes messengers, representing the photons, 
through a network of message-processing units, playing 
the role of the optical components in the laboratory ex- 
periment. In what follows we give a detailed description 
of each of the components of the network representing 
the complete experimental setup of the quantum eraser 
experiment, schematically depicted in Fig. [T] 



A. Messenger 

A messenger has its own internal clock, the hand of 
which rotates with frequency /. When the messenger is 
created, the hand of the clock is set to time zero. As the 
messenger travels from one position in space to another, 
the clock encodes the time-of-flight modulo the period 
1/f- The message, the position of the clock's hand, is 
most conveniently represented by a two-dimensional unit 
vector e ; = (e j,eij) = (cos ipi , sin ipi ) , where ipi = 2nft, 
the subscript I > labeling the successive messages. The 
messenger travels with a speed c/n where n is the re- 
fractive index of the medium in which the messenger 
moves and c is the light velocity. Clearly, this messen- 
ger is the event-based equivalent of a classical, linearly 
polarized electromagnetic wave with frequency /: The 
messenger corresponds to the light ray with wave vector 
k(fc = 2%f/c) and the clock mimics one of the electric 



B. Beam splitter 

Here we construct a processing unit that acts as a BS, 
not by calculating the amplitudes according to QT, but 
by processing individual events (see Fig. [2]). It consists of 
an input stage, a simple deterministic learning machine 
(DLM)ii±&±&22i22, a transformation stage (T), an out- 
put stage (O) and has two input and two output channels 
labeled by k = 0, 1. We now define the operation of each 
stage explicitly. 

• Input stage: The DLM receives a message on ei- 
ther input channel or 1, never on both chan- 
nels simultaneously. The arrival of a message on 
channel (1) is named a (1) event. The input 
events are represented by the vectors v; = (1,0) or 
v; = (0, 1) if the Zth event occurred on channel 
or 1, respectively. The DLM has six internal regis- 
ters = {C ktl ,S^ t ), Y^ ; = (C^^S^j), = 
(Cjf;, S^i) and one internal vector x; = (xq,i, 
where xo,; + x i,i = 1 an d x^j > for k = 0, 1 
and ahW > 0. These seven two-dimensional vectors 
are labeled by the message number I to indicate 
that their values may change every time the DLM 
receives a message. The DLM has storage for no 
more than fourteen numbers. 

Upon receiving the Ith input event, the DLM per- 
forms the following steps: It stores the first two 
elements of message ykj in its internal register 
^ki = i^kl' ^hi), tne m iddle two elements of yk,i 
in Y^ ; = (C^^S^t), and the last two elements 
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of yk,i in = (C£ l: SjT/)- Then, it updates its 
internal vector according to the rule 

x; = 7x ; _i + (1 - 7)vj, (8) 

where < 7 < 1. Note that by construction 
^o,i + x i,i = 1, 2:0,; > and > 0, and the 
DLM stores information about the last message 
only. The information carried by earlier messages 
is overwritten by updating the internal registers. 



From the solution of Eq. (JH) , 

1-1 

x, =7 i x ; _ 1 + (l- 7 )]Ty--'- 1 v J+1 , (9) 
j=o 

the fact that in practice the sequence 
{vi, V2, • • • ,vk} is finite, and the usual trick 
to assume a periodic continuation of the sequence, 
we have 



1K-1 



K-l 



X,„K ^ l K * (m -l)K + (1 ~ 7) 3 S+l = l K *(rn-l)K + (1 - 7) ^ ' lv i+l+(m-i; 

j=( m -l)K j=0 

K. 



K 



7 *-(m-i)K + (1 - 7) f x 



r 



(10) 



where 

and m > 0. From Eq. (|10|) we find 

*mK = 7 mK X + (1 - 7) . 7 „ f, 



forms them into a new eight-dimensional vector 



1-7 



A" 



K, 



and hence 



lim x mK - = 2 

m^oo 1 — 7^ * — ' 



i=o 



such that 
lim 



1 A ' _1 

n lim Xmif = — v J+ i 



i=0 



(11) 



(12) 



(13) 



(14) 



From Eq. (fT4|) . we conclude that as 7 — * 1~ the 
internal vector converges to the average of the vec- 
tors Vi , V2 , • • • , vk which represents the relative 
frequency of input events at the two channels of 
the BS (k = 0,1). The parameter 7 controls the 
speed of learning and also limits the precision with 
which the internal vector can represent a sequence 
of constant input messages^!. Disregarding the fact 
that according to Eq. (|14p . we should let 7 — > 1~ 
to obtain the limiting value of the average of the 
v's, it is the only free parameter in the model. In 
practice, in the simulation we fix it once and for all. 

• Transformation stage: The second stage (T) ac- 
cepts the messages from the input stage, and trans- 



(15) 



1,1^1,. 
oh r^p 

CV qP 

qV qP 
°0.l D 0.l 
qH r<P 

qH r<P 

CV qP 
°0.l J 0,l 
qV Ql 



If we rewrite the transformation T using complex 
numbers, we find that 





bf 




f a o\ 

v 

«0 

a? 

v«ry 



(16) 



which is the unitary transformation in the quantum 
theoretical description of a BS, if (a^ ,Oq ,aY) 
and (&o , by, bf, b\) denote the input and output 
amplitudes of the photons with polarization H and 
V in the and 1 channels of a BS, respectively. 
Note that in our simulation model there is no need 
to introduce the (quantum theoretical) concept of a 
vacuum field, a requirement in the quantum optical 
description of a BS. 

Output stage: The final stage (O) sends out a mes- 
senger (representing a photon) carrying the mes- 
sage 



wi,i/s ,i 

W2,l/Slj 

W3,l/Sl,l 
SQ,l/s 2 ,l 

V si.i/s 2 ,i I 



(17) 
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where 
w .i 

11)2,1 

SQ,l 
Sl,l 
S2I 



n H n p p? — qH r <p r- — 

u cy°(yv a ' > i — D i,i^i,i\ x i-h 

n v q p i q v qP i 

^o,i^o,iV Xo ^ ~ °i,l°i,i\ x i,l> 

riV qP i i qV qP i 

{ ^i,i D i,iV Xl ' 1 ' °o,i a o,iv x o,h 



J 0,l 



J 1,V 



J 2,l 



w 



3,1' 



J 0,l 



J ll 



"2,1 



u 3,l' 



(18) 



through output channel if s\ l > 2r where < r < 
1 is a uniform pseudo-random number. Otherwise, 
if s\ l < 2r, the output stage sends through output 
channel 1 the message 



where 

zo.i 

Zl.l 
z 2.l 

Z3,l 

to,i 



( zo,i/to,i \ 
z 2,i/h,i 

Z3,l/tl,l 
to,l/t2,l 

V h,i/t 2 ,i ) 



r H r p p? — oH n p r- — 

^i,i u i,W x i,i ~~ * 5 cu o,zv a;o . i ' 

^o,fio,lV x o,i + Si,fii,i\f x ~ij, 

n v q p i q v q p i 

°l,i D l,iV a ' 1 ) i — ^cm^o.zv x o,h 



(19) 



zl , + z 2 



■o,i 



1,1' 



■2,1 



-3,1' 



■0,1 



2.1 



-3,1- 



(20) 



The use of pseudo-random numbers to select the out- 
put channel is not essential^. We use pseudo-random 
numbers to mimic the apparent unpredictability of the 
experimental data only. Instead of a uniform pseudo- 
random number generator, any algorithm that selects 
the output channel in a systematic manner might be em- 
ployed as well™. This will change the order in which 
messages are being processed but the content of the mes- 
sages will be left intact and the resulting averages do not 
change significantly. 



C. Polarizing Beam Splitter 

A polarizing beam splitter (PBS) is used to redirect 
the photons on the basis of their polarization (H or V). 



The structure of the event-based processor that simulates 
a PBS is identical to the one of the BS and differs in the 
details of the transformation stage only. For the PBS, 
the transformation T reads^i 



°0,Z + 1 
oH 

D 0,l+1 

OA 

J l,l+1 

°o,i+ 
r v 



C ([l+W x 0,l+1 
C 0,l+ly/ x 0,l+l 

i s i,i+ W x u+i 
s i,i+i\/ x id+i 
c u+W x u+i 



\ 



c u+W x u+i 

i^o,l+iy/ x °' l+1 
Sq,i+W x o,i+i / 



(21) 



D. Remaining optical components 

In contrast to the BS and PBS, in terms of message 
processing the HWP and QWP are passive devices in 
the sense that the adaptive unit, the DLM, is not re- 
quired for a proper functioning of the devices. As can be 
seen from the quantum theoretical description (see Ap- 
pendix), a HWP does not only change the polarization 
of the photon but also changes its phase and a QWP 
additionally, introduces a phase difference between the 
H and V components. In our simulation model, the 
functionality of these optical components is implemented 
through plane rotations of the vectors (cos^;, sin^;), 
(cos ipgj, sin and (cos ip%j, sin ip^). 



E. Data gathering and analysis procedure 

In the simulation, the data is collected in the same 
manner as in the experiment. Detector Do (Di) registers 
the output events at channel (1) (see Fig. [TJ. During a 
run of TV events, the algorithm generates the data set 



r = {x l \l = l,...,N;4r,6 



HWPO > °H WPl j "QWP 



}, (22) 



where xi = 0, 1 indicates which detector fired (D or Di), 
<j) denotes the phase shift (proportional to the difference 
in time-of-flight of PathO and Pathl) between the two 
interferometer arms and Qhwpo, Ohwpi, 9qwp denote 
the angles of the optical axis of the respective waveplates 
with the laboratory frame. For fixed Ohwpo, Qhwpii 
^ QWP and <t>, the number of detection events in detector 

1 is given by N\ = Yli=i x i anc ^ = N — N\ is the 
number of detection events in detector 0. The appearance 
of interference fringes is conveniently characterized by the 
visibility^ 



V = 



N max N n 



(23) 



where N max and N m i n denote the maximum and mini- 
mum of N for all <j) £ [0, 2tt[. Notice that for the ex- 
periment depicted in Fig. [TJ the visibility is a function of 
&hwpo, Qhwpii and 0qwp- 
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VI. SIMULATION RESULTS 

The processing units that simulate the optical com- 
ponents are connected in such a way that the simula- 
tion setup is an exact one-to-one copy of the real experi- 
ment (see Fig.[T]). The simulation procedure is as follows: 
For each choice of <p in the range [0, 2ir[, we fix Ohwpo, 
Ohwpi and Oqwp and perform a simulation with 10 6 
events, randomly distributed over groups of Nh = 200 
or Ny = 200 events («i = H and = V in the notation 
of Section HVj). Then for each choice of Ohwpo, Ohwpi, 
Oqwp, we repeat this procedure. The result of these cal- 
culations form the data set T (see Eq. (|22[) ). From this 
data set, we compute the visibility according to Eq. (f23|) . 
All simulations have been carried out with 7 = 0.99. 



A. Without QWP 

In Fig. [3] we show our simulation results for the vis- 
ibility as a function of 20hwpi for the quantum eraser 
experiment with the QWP removed (see Fig. [T]). First 
we consider the case in which the source emits photons 
that in QT are described by a pure, vertically polarized 
(V) state. Each such photon, after passing through the 
first BS, has equal chance to end up in either of the 
two arms of the interferometer. In our simulation, the 
messenger representing this photon carries the message 
(O,O,cos^,sin'0o r ,O, l) (see Eq. ©). If the photon 
follows PathO, it encounters HWP0, the optical axis of 
which makes an angle Ohwpo with respect to the labora- 
tory frame. HWP0 rotates the polarization of the photon 
by an angle 20hwpot^- The event- by-event simulation 
data and the results of QT are shown in Fig. EJa). The 
simulation data are in quantitative agreement with the 
averages calculated from QT and in qualitative agree- 
ment with the experimental data (see Fig. 4(a) in Ref4). 

Next, we consider the case where in QT, the input to 
the quantum eraser is described by a (completely) mixed 
state. In QT, a mixed state simply means that photons 
emitted by the source are described by an incoherent mix- 
ture of horizontally and vertically polarized pure states. 
In Section IIV1 we explained how to implement mixed 
states in the event-based simulation approach. The sim- 
ulation data for a source emitting photons described by a 
(completely) mixed state are shown in Fig.^b) and (c). 
Also in this case, our simulation data are in quantitative 
agreement with the averages computed from QT and in 
qualitative agreement with the experimental results re- 
ported in Ref£ (see Fig. 4(b) and (c)). 



B. With QWP 

In Fig. H]we present some simulation data for the case 
that the QWP is present, see Fig. [TJ and 9qwp — 0. If 
Oqwp = 0, the QWP does not change the polarization 



1 ■ ■ 

_ w / 


I - r 1 1 . ■ T • _ i ■ ■ • ■ 

\k \ s ■ \ 
\ ■ ; > / t \ 




\ k r y ? \ - 






(iii 


1 1 V 1 k 1 1 1 1 1 



15 30 45 60 75 90 105 120 135 150 165 180 




15 30 45 60 75 90 105 120 135 150 165 180 




15 30 45 60 75 90 105 120 135 150 165 180 



FIG. 3: Visibility as a function of the angle 29hwpi for 
the quantum eraser experiment with the QWP removed (see 
Fig. [T}. The markers (squares, bullets) and lines (solid, 
dashed) represent the event-by-event simulation data and the 
quantum theoretical results (see Eqs. (|30|) and (|3ip ). respec- 
tively, (a) The source emits photons described by the pure 
vertically polarized state V and 9hwpo = 45° (red bullets and 
solid line), Ohwpo = 10° (black squares and dashed line); 
(b) The source emits photons described by the completely 
mixed state (pv = ph = 1/2) and Ohwpo = 45°; (c) The 
source emits photons described by a partially mixed state with 
pv = 2/3, ph = 1/3 and Ohwpo = 22.5° (black bullets and 
solid line). The red dashed and blue dotted curves represent 
the quantum theoretical results for the pure vertically polar- 
ized state V and the completely mixed state, respectively. 



of the photons but changes their phase. We only con- 
sider the case that the single-photon source emits pho- 
tons that in QT, are described by a pure state. Fig- 
ure [UJa) shows simulation data corresponding to incom- 
ing V-polarized photons, for Ohwpo — 45° (red bul- 
lets) and Ohwpo = 10° (black squares). In Fig. rjjb) 
we show the simulation data for the source that emits 
photons in a state that QT would characterize with 
£ = 45°, and Ohwpo = 22.5°. In our simulation, this 
state is represented by messengers that carry the mes- 
sage (cos tpQ , sin ipQ , cos 4>o > sm > 1 / v2, 1 /V%) ■ As in 
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FIG. 4: Visibility as a function of the angle 26hwpi for the 
quantum eraser experiment depicted in Fig.[T]with Qqwp = 0. 
The markers (squares, bullets) and lines (solid, dashed) rep- 
resent the event-by-event simulation data and the quantum 
theoretical results (see Eq. (|32[) ). respectively, (a) The source 
emits photons described by the pure vertically polarized state 
V and Ohwpo = 45° (red bullets and solid line), Ohwpo = 10° 
(black squares and dashed line); (b) The source emits pho- 
tons described by the pure £ = 45°-polarized state and 
Ohwpo = 22.5°. 



all other cases shown, the agreement between the event- 
based simulation data and QT is excellent. 



tons in our simulation model have full which-way infor- 
mation. A consequence of our model is thus that the 
relation V 2 + T> 2 < 1 cannot be regarded as quantifying 
the notion of complementarity: Our model always allows 
a particle-only description of the quantum eraser exper- 
iment, independent of the purity of the state describing 
the photons in quantum theory. 

In summary, concepts of quantum theory applied to 
individual events fail to provide a logically consistent ex- 
planation for the experimental observation of single de- 
tector "clicks" building up an interference pattern and 
leave no option but to postulate that "this is the way 
it is". In contrast, our event-based simulation model, a 
classical locally causal dynamical system, reproduces the 
results of quantum theory without making reference to 
the solution of a wave equation and provides a simple, 
particle-based mental picture for what each individual 
photon experiences as it travels from the source to the de- 
tector. Just like in the experiments, our model produces 
data sets Eq. (|22|) which can be given to a third party 
for analysis long after the simulation has been finished. 
Because of the strong similarity between the experimen- 
tal and simulation data sets the third party will have a 
very hard time, if possible at all, to identify the data sets 
as originating from a so-called "quantum experiment" or 
from a "classical simulation model" . 

Finally, we would like to emphasize that the algorithms 
used to simulate the optical components of the quan- 
tum eraser have not designed to exclusively simulate this 
particular example but they can be used to reproduce 
the results of many other quantum optics experiments as 

„, C Y [ 17. 18. 19. 20. 21. 22. 23. 24.25. 26. 27, 28. 29 



VII. DISCUSSION 

We have demonstrated that our classical, locally 
causal, particle-like simulation approach reproduces the 
results of the quantum eraser experiment^ and the re- 
sults of quantum theory describing the averages of these 
experimental results. 

During the event-by-event simulation of the quantum 
eraser experiment we always have full which-way infor- 
mation of the photons (messengers) since we can always 
track them. Nevertheless, depending on the settings of 
the optical apparatuses, the photons build up an inter- 
ference pattern at the detector. Although the appear- 
ance of an interference pattern is commonly considered 
to be characteristic for a wave, we have demonstrated 
that, as in experiment, it can also be built up by many 
photons. These photons have full which-way informa- 
tion and arrive one-by-one at a detector. Hence, even in 
the case that the source emits single photons, described 
by a pure state in quantum theory, and that V = 1, 
commonly associated with full wave character, the pho- 



APPENDIX 

According to quantum theory (QT), photons in a pure 
state are described by the state vector 



f a o\ 



V 

1 



(24) 



where H and V refer to the horizontal and vertical direc- 
tion of polarization and the subscripts refer to the wave 
in PathO and Pathl, respectively. Within QT, the ac- 
tion of the various optical components is defined by the 
matrices 



Tbs = —7= 




(25) 
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T QWP (6) = — | „ „ , ; | . (20, 

Tpss - I n n i n l> (26) 
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Thwpo{0) = -i \ n n i n , (27) 



T, 



HWPl 



(0) 



10 
10 

c s 
s -c 



where 9 denotes the angle of the optical axis with respect 
(28) to the laboratory frame, c = cos 26 and s — sin 26. 

I 



Using these expressions, it is somewhat tedious but straightforward to calculate the visibility Eq. ([2"5)> . We list the 
expressions for the cases for which we perform event-based simulations. 

1. With the QWP removed, see Fig.[U and for incoming photons that are described by a pure state of polarization 



V 



where 6q = 6hwpo and 6\ = 6hwpi- 



2 sin(£ - 26 + 26>i) sin(£ - 26>i) 



sin 2 (£ - 26»o + 26»i ) + sin 2 (£ - 26>i) 



(30) 



2. With the QWP removed and for incoming photons described by a mixed-state photon input with pv/ph = 
tan 2 /3: 



V 



2sin(26» - 26»i)sin(2(9i) + 2 tan 2 /3cos(26» - 20i) cos(20i) 



sin^(26» - 20i ) + snT(20i) + tan 2 /3[cos 2 (2(9 - 26>i) + cos 2 (26»i)] 



(31) 



3. With the QWP present, 6qwp = 0°, and for incoming photons that are described by a pure state of polarization 



V 



[sin 2 46»i sin 2 (2£ - 26 ) + 4[cos 2 26>i sin(£ - 26 ) sin£ - sin 2 26 x cos(£ - 26 ) cos^] 2 ] 



1/2 



cos 2 26»i sin 2 (£ - 29 ) + sin 2 26>i cos 2 (C - 26 ) + sin 2 26»i cos 2 ^ + cos 2 26>i sin 2 ^ 



(32) 
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